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Abstract 

The structure factor for a Bose system with Bose-Einstem condensate is consid- 
ered in the frame of the self-consistent mean-field approximation. The accomplished 
analysis demonstrates the principal importance of the following three points: the 
necessity of preserving the approximation order, the necessity of taking into ac- 
count anomalous averages, and the necessity of gauge symmetry breaking. If any 
one of these necessary conditions is not satisfied, calculations yield the appearance 
of unphysical divergences of the structure factor, which implies the occurrence of 
fictitious instability, which contradicts experiments. 
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1 Introduction 



The structure factor is an important quantity characterizing, as follows from its name, 
the structural properties of the considered system as well as the spectrum of collective 
excitations. The structure factor can be measured in scattering experiments. It is also 
related to other directly measurable quantities, such as the isothermal compressibility and 
sound velocity. 

The structure factor of Bose systems has been repeatedly investigated for the case of 
liquid ^He, both theoretically [1,2] and experimentally [3-9]. In recent years, the physics of 
dilute Bose gases [10] has attracted a great deal of attention (see review articles [11-15]). 
The structure factor of cold trapped atoms was also measured [16]. Contrary to liquid 
helium, for which strong interparticle correlations must be accurately taken into account 
[2,17], dilute gases allow for their description the use of the mean-field approach. Thus, 
the measured structure factor of dilute trapped ^^Rb atoms was found [16] to be in good 
agreement with the Bogolubov approximation [18,19]. The latter is known to be well justi- 
fied for weakly nonideal Bose gases, with asymptotically small interactions. When atomic 
interactions are not asymptotically weak, but rather are finite, their theory encounters the 
notorious Hohenberg-Martin [20] dilemma of conserving versus gapless approaches. This 
dilemma has recently been resolved by constructing a representative statistical ensemble 
for Bose systems with broken gauge symmetry [21,22] and by developing a self-consistent 
mean-field approximation [23-26]. 

The aim of the present paper is to consider the structure factor of a dilute Bose- 
condensed system, in the frame of the self-consistent mean-field approximation [21-26], 
and to analyze the properties of the structure factor. This approximation, as is shown in 
Ref . [26] , can be employed for both weak as well as strong atomic interactions. Under the 
Bose-condenscd system, one can keep in mind atomic dilute Bose gases [10-15], molecular 
systems formed by bound bosons [27-30] and fermions [31-35], which can be achieved by 
invoking the Feshbach resonance [36,37], or other effectively bosonic assemblies. 

Throughout the paper the system of units is used, where the Planck constant and 
Boltzmann constant are set to unity, h = 1, kB = 1. 



Phase transitions are known to usually be accompanied by spontaneous symmetry break- 
ing [38]. The appearance of Bose-Einstein condensate is associated with the breaking of 
global gauge symmetry. The most convenient method of breaking the latter is by means 
of the Bogolubov shift [39,40]. Then, the field operator is represented in the shifted form 



in which i]{r,t) is a nonoperator function, defining the condensate wave function, and 
■(/'i(r, t) is the field operator of uncondcnscd particles, satisfying the Bose commutation 
relations. The field variables of condensed and uncondensed particles are orthogonal to 
each other. 



2 Bose-condensed system 





(2) 
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The condensate function is normalized to the number of condensed particles 



No = l \vir,t)\Ur, (3) 
while the operator of uncondensed particles satisfies the normalization 

TVi ^<Ni>, (4) 
given by statistically averaging the operator 

7Vi = I i;l{r,t)M^,t) dr (5) 

for the number of uncondensed particles. The average number of particles in the whole 
system is N — Nq + Ni. The statistical condition 

<^i(r,i) (6) 

guarantees the conservation of quantum numbers. 

Gauge symmetry can also be broken by introducing infinitesimal sources and defining 
the Bogolubov quasiaverages [39,40]. Both these methods of gauge symmetry breaking, 
by means of the Bogolubov shift (1) and by means of infinitesimal sources, have been 
shown to be equivalent and sufficient for describing Bose- Einstein condensation [39-45]. 

The operators of observables, after the Bogolubov shift (1), are defined on the Fock 
space generated by ^| (see mathematical details in Refs. [23,46]). The energy 

operator is taken in the standard form 

H = j V^^(r) ^ + V'(r) ^^+\ j i^K'^Wi'^')^^^ - rOV'lrOV^lr) dr dr' , (7) 

in which, for breavity, the time dependence of the field operators is not shown explicitly; 
U — U{r,t) is an external field, if any; and $(r) is the interaction potential. 

In what follows, we shall consider a uniform system, with no external fields, U — 0. 
The system is assumed to be equilibrium. An equilibrium statistical system is to be char- 
acterized by a representative statistical ensemble [47], which is the couple {^{1^1), p} of 
the Fock space J-'{ipi) and the appropriate statistical operator p. The latter can be found 
from the minimization of the information functional [46,47], incorporating the normaliza- 
tion conditions (3) and (4). This yields [21-26] the statistical operator 

(8) 



^ Tre-/^^['?''/'il ' 

in which (5 — 1/T is inverse temperature and the grand Hamiltonian is 

H[r],^/ji]^H-l^oNo-l^iN,, (9) 

where /iq and pi arc the Lagrange multipliers guaranteeing the validity of the normaliza- 
tion conditions (3) and (4). 
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For a uniform system, we may set 

with V being the system volume. From $(r) and we pass to their Fourier transforms 

$k and Ok, respectively. The lowest nonzero averages of the products of are the normal 
average 

nk = < OkOk > (11) 

and the anomalous average 

(7k = < akO-k > ■ (12) 

The average < ak >= is zero, according to condition (6). 

The higher orders of the products of ak can be simplified in the Hartree-Fock-Bogolubov 
approximation [22-26]. To this end, we define 

UJk = — + P% + Po$k + 77 ^P^k+p - 1^1 (13) 

and 

Ak = Po$k + E '^P^k+p . (14) 

^ p^O 

Then the grand Hamiltonian (9) acquires the form 

Hhfb = Ehfb + E <^kakak + ^ E ^k {aW-k + Q-kOk) , (15) 
with the nonoperator part 

Ehfb = Po$o - l^o) Nq - I- pl<^>oV - ^ E ^k+p(nknp + cTkCTp) , (16) 

where pi = Ni/V. 

The grand thermodynamic potential is 

Q = -rinTre-^-^t"''^!! . (17) 

For thermodynamic stability, the number of condensed particles Nq is to be a minimizer 
of the grand potential (17), which implies 

The latter defines the condensate chemical potential 

/^O = P$0 + E (^P + ^P)*P ■ (19) 
^ p^O 
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The grand Hamiltonian (15) can be diagonalized by means of the Bogolubov canonical 
transformation 



Ok = Mk&k + ^^-k^-k , % = , % = , (20) 

2£k 2£k 

with the Bogolubov spectrum 

£k = . (21) 

Then Hamiltonian (15) reduces to the diagonal form 

Hb^Eb + Y1 ^Aby. , (22) 

in which 

Eb = Ehfb + - E (^k - c^k) • (23) 

^ k^O 

The existence of Bose- Einstein condensate requires the condition 

lim£k = 0, £k>0, (24) 
fe— >o 

from which the chemical potential of uncondensed particles follows as 

/xi = p$o + E (^p - ^p)^p ■ (25) 

The Lagrange multipliers (19) and (25) are, evidently, different. They would coincide 
only for the ideal gas, when $k = 0, hence /xq = A*i = 0. But for a nonideal system. 
Ho Hi- Sometimes, one sets these potentials equal by omitting the anomalous average 
CTp, which corresponds to the Shohno model [48]. However, direct calculations [25,26,49] 
demonstrate that the anomalous average is of the order and can even be larger than the 
normal average. Therefore, omitting the larger quantity than what is retained has nothing 
to do with a reasonable approximation. In addition, neglecting the anomalous average 
makes the system unstable [21,50,51]. 



3 Structure factor 

The structure factor is an important characteristic, which can be measured in scattering 
experiments, such as neutron or photon scatterings [52]. Let us consider, for concreteness, 
neutron scattering. A neutron, with the initial energy Ei and momentum kj, scatters on 

an atom, acquiring the final energy Ef and momentum k/. The momentum and energy, 
transferred by a neutron to an atom, are k = kj — kj and uj = Ei — Ef. The double- 
differential scattering cross section of neutrons on the system can be written [38] as 

1 (Pa o kf ^ , 

iVdii;SJ = ''2*rS«(k.-). (26) 
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where fik is a solid angle around k; bs is the scattering length of the neutron on an atom; 
and StotO^, is the total dynamic structure factor. The latter is defined by the expression 

StotO^, ~ ]y J -^(^' ^' ^'^ ^) ^-^P {~'^k(r — r') + iujt} drdv'dt , (27) 

in which 

R{y, t, y', f) = < n(r, t)h{r\ t') > (28) 
is the density-density correlation function, and the density operator is 

n{r,t) = i^\r,t)i^{r,t). (29) 

The total dynamic structure factor (27) can be decomposed into two terms 

StotiK uj) = 'S'e(k, u) + S{k, u) ; (30) 

the first term 

Se{k, uj) = ^ J p{r, t)p{r', 0) exp {-ik(r - r') + iuut} drdr'dt , (31) 

with p(r, t) =< n(r, t) >, is responsible for elastic scattering and is called the disconnected 
dynamic structure factor; the second term, called connected, 

S{k, uj) = ^ J [R{t, t, r', 0) - p(r, t)p{r', 0)] exp {-ik{r - r') + iujt} drdr'dt (32) 

describes inelastic scattering. 

The elastic dynamic factor (31) is nontrivial for nonuniform systems [53,54], but for a 
uniform equilibrium system, when p(r, t) = p, it is very simple, being 

,Se(k,u;) = (27r)V(k) %) . (33) 

Therefore, for uniform systems, more attention is payed to the inelastic dynamic factor 
(32). The latter enjoys the sum rules, one of which 

I-oo '^(k>^)^ = '^(k) (34) 
defines the static structure factor 5'(k), another gives 

and the third, connects the integral 

/■~ 1 du; Rex(k, 0) 

/ - S(k, u;) — ^ — , 36 

J-oo uj 27r Ip 

with the responce function x(k, a;). At very large momentum k — > oo, there are the 
following asymptotic equalities: 

du) 



/oo (ii,] 
_^ 5(k,a;)— :^1 (A; ^ oo) , 
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/ \U— ] S[K,CaJ) — ~0, 

r ( V 2e <k> 



27r 3m N 



here < > is the average kinetic energy of the system. 
For an equihbrium system, where 

p{r,t) = p{r) =<n(r)>, (37) 

with h{r) = n(r, 0), the static structure factor, defined in Eq. (34), becomes 

S{k) = -i y [< n(r)n(r') > - p(r)p(r')] exp{-iA;(r - r')} drdr' . (38) 

Let us introduce the Fourier transform of the density operator (29), 

Pk = / n{v) e-"^-^ dv , n{v) = :^ E Pk e^'^'^ ■ (39) 

The following properties hold: 

[Pk, Pp] = 0, [n(r), n(r')] = , p,^ = p_k , n+(r) = n(r) . 
The zero component of pk yields the total number of particles, 

po = y" ri(r) dr = N . (40) 

The static structure factor (38) can be expressed as 

= ^ (41) 

through the dispersion 

A'(pk) = < pl^Pk > - < p^ >< pk > . (42) 
The central value of the structure factor (41), 



S{0) = = pTKT = (43) 



connects it with the particle dispersion A^(7V) = A^(po), the isothermal compressibility 



where P is pressure, and with the hydrodynamics sound velocity s, 

(45) 

Thus, the isothermal compressibility and sound velocity are expressed through the struc- 
ture factor as 

Relations (43) and (46) are exact, being valid for any equilibrium, not necessary uniform, 
system. 

Passing to a uniform system, for the Bogolubov shift (1), we have 

V'(r) ^Vp~o + V'i(r) . (47) 
The density operator (29) becomes 

n(r) = po + ni(r) + Vp^[V'I(r)+V'i(r)] , (48) 
where the density operator of uncondensed particles is 

ni(r) = V'I(r)V'i(r) . (49) 
The Fourier transform of the total density operator (48) takes the form 

Because of the properties 

< aj^Op > = SkpUk , < Ok > = , 

we have 

< Pk > = SkoN . (51) 

For dispersion (42), we find 

A2(pk) = AT + 2No{nk + (Jk) + ^ < a^aj,+qak+paq > - ^oA^i' ■ (52) 

p,q7^o 

When the interaction potential $(r) is singular at r ^ 0, then pair atomic correlations 
are very important [52,55]. But here, we have assumed that the interaction potential 
is soft, so that it possesses the Fourier transform $k- For the soft potential, we can 
apply the Hartree-Fock-Bogolubov (HFB) approximation. Employing the latter, one has 
to be cautious, since its blind application can lead to unphysical results. Thus, if we 
would formally decouple the four-operator term in dispersion (52), by using the HFB 
approximation, we would come to 

A'(pk) — >N + 2Noin^ + a^) + ^(nk+pnp + (7k+p(7p) . (53) 

Pt^O 
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However, one should keep in mind that the use of the HFB approximation results in 
the Hamiltonian (22), which is quadratic with respect to the operators of uncondensed 
particles. Therefore, all expressions of the order higher than two with respect to these 
operators are not well defined in this approximation. One could consider such high-order 
terms only if their behavior were reasonable. But sometimes these higher-order terms 
result in unphysical divergencies, which simply means that they should be omitted, being 
outside of the region of applicability of the used approximation. In the second-order ap- 
proximation, only the second-order terms are well defined [14,21,23,50,51]. Consequently, 
dispersion (52), in the HFB approximation, has to be written as 

A2(pk) = A^ + 27V(nk + ak). (54) 
Recall that the normal and anomalous averages in the HFB approximation [23-26] are 

1 „ , , ^^^^ 



rik = coth 

2£k 



— - , CTk = coth 

2TJ 2 ' 2£k 



2T 



where 6(;k, Ak, and £k are given by Eqs. (13), (14), and (21), respectively. The Bogolubov 
spectrum (21), according to condition (24), is gapless, having the phonon form £k — ck 
at k ^ 0, with the sound velocity c defined by the equation 



T7 E ^P^P • 



Dispersion (54) is finite. For instance, in the case of k 
we have 

A^(7V) = N 



(56) 

0, when A2(po) = A'^{N), 



1 + 2 lim(nk + (7k) 



Using Eqs. (55), we get 



Therefore, 



lim(nk + CTk) 

K— »0 



T 



- 1 



NT 



Consequently, 



5(0) = 



T 



r.2 ' 



s — c . 



The structure factor (41), according to Eq. (54), takes the form 

>S(k) = 1 + 2(nk + CTk) . 
Substituting here expressions (55), we find 

fc2 



5(k) 

Note that the Feynman relation 



coth ( ^ 

2m£k V2T 



(57) 
(58) 

(59) 
(60) 



2m5(k) 



(r = o) 
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is valid only at zero temperature. The asymptotic properties of the structure factor (60) 
are: 



S{\i) ~ 1 (A; ^ oo) . 

If instead of expression (54), one would use Eq. (53), then the structure factor 5'(k) 
would be divergent, which has no physical sense. Hence, Eq. (54) provides the correct 
mean-field approximation. The correctness of Eq. (54) will be analyzed in more detail in 
the following section. 



4 Approximation analysis 

The form of the structure factor (60) is due to that of dispersion (54). Here, we shall 
analyze the correctness of the approximation in which Eq. (53) has been replaced by Eq. 
(54). 

We shall also stress the necessity of correctly taking into account the anomalous aver- 
ages and the importance of gauge symmetry breaking. 



4.1 Preservation of approximation order 

The replacement of Eq. (53) by Eq. (54) is based on the requirement of preserving the 
approximation order with respect to the operators of uncondensed particles. Suppose 
that we would go outside of the region of applicability of the considered second-order 
approximation and would retain the last term in Eq. (53). Then A^(pk) would be 
proportional to iV^/^, hence <S'(k) proportional to A?"^/^, that is, the structural factor S'(k) 
would be divergent, hence, unphysical. To show this, consider A^(po) = A^(A^). Then 
the last terms in Eq. (53) give 

J (27r)3 27r2 ' J (27r)3 277^ 

Therefore, 

TT 

In this way, when T > 0, then in thermodynamic limit 

S{fS) oc N^/^ ^ oo , s oc A^-^/^ ^ , 

N^l^ ^ oo (TV 

This means an absolute instability of the system at any finite temperature. 
At zero temperature, Eqs. (55) reduce to 

l + 2g^ 1 1 f _ k 
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Prom here, 

lim(nk + = - ^ (7^ = 0). 

The last terms in Eq. (53), at T = 0, are finite, since 

[ o dli f n dk 

The dispersion A^(7V) (x N is thermodynamically normal. The value S{0) oc const > is 
finite. However, as follows from relations (46), 

1 

ktoc — — >-oo, socT— >-0 (T— >-0), 

which tells that the system is again unstable. Thus, at any temperature T > 0, in 
thermodynamic limit, one has s = and kt — oo. Hence, such a system would always be 
unstable. But, of course, the appearance of this instability is not physical and is caused 
only by going outside of the region of applicability of the used approximation. Preserving 
the approximation order requires to reject the last terms in Eq. (53). Then no unphysical 
divergencies appear. 



4.2 Importance of anomalous averages 

One of the popular tricks when considering Bose-condensed systems is the omission of the 
anomalous average (Jk, which corresponds to the Shohno model [48]. One often ascribes 
this unjustified trick to Popov. But, as easy to check from the original works by Popov 
[56,57], which are cited in this respect, he has never suggested this trick. Here we show 
that omitting the anomalous average is principally wrong, since then the system becomes 
unstable and its structural factor infinite, which is, certainly, of no sense. 

Assume that the anomalous average Uk in Eq. (54) is set to zero. Then we have 



At nonzero temperature. 



Since k ~ we get 



A2(iV) ^ N (l + 2\iuink^ 



TflT 

^^--^ (T > 0, A; ^ 0) . 



A\N) ^ ^ N'/' (T > 0) . 

pZ/6 



Then, at finite temperature. 

These quantities diverge in thermodynamic limit, which implies instability. 
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At zero temperature, 

TTIC 

rik ~ 2^ (T = 0, A; ^ 0) . 

Prom here 



Hence, as T — > 0, we have 



5(0) ~ ^ iV'/3 , . . 

The structural factor diverges in thermodynamic limit, and the isothermal compressibility 
Kt becomes divergent even at finite N, when T — > 0. 

Thus, at all temperatures T > 0, we come to the senseless unphysical behavior, when 

5'(k) ^oo, s ^ , Kt —>■ oo . 

This proves that the anomalous average plays a principally important role and cannot be 
omitted. 



4.3 Gauge symmetry breaking 

Gauge symmetry breaking is known to be sufficient for describing Bose-Einstein con- 
densation. This follows from the general theorems, independent from approximations, 
proved by Bogolubov [39,40], Ginibrc [41], and Lieb et al. [42-45]. At the same time, 
in the Penrose-Onsager [58] scheme of Bosc condensation no explicit breaking of symme- 
try is required, but the appearance of condensate is manifested in the properties of the 
density-matrices eigenvalues or their order indices [59,60]. Here wc give arguments that, 
in the frame of a mean-field approximation, gauge symmetry breaking is necessary. 

Suppose, wc consider a mean-field approximation without breaking the gauge symme- 
try. Then the numbcr-of-particlc dispersion has the standard form 

A\N) = ^nk(l + rik) = A\No) + A\N^) , 

k 

which can be separated into two terms, 

A2(7Vo) = No{l + No) , A2(7Vi) = ^ n^{l + n^) , 

corresponding to condensed and uncondensed particles, respectively. The occurrence of 
Bose-Einstein condensation implies that the number of condensed particles becomes pro- 
portional to the total number of particles, Nq (x N. As a result, 

A^(No) oc N^ (No oc N) . 

Consequently, A^(iV) oc N^. Therefore, the isothermal compressibihty and the structure 
factor diverge in the thermodynamic limit as 

5'(k) ocA^— >-oo, Kt ex. N ^ oo , 
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which means instabihty. 

In this way, the use of the mean-field approximation requires that gauge symmetry 
be broken. It may be that the exact consideration or some refined approximations (see, 
e.g., Ref. [17]) allow for a description without gauge symmetry breaking. However, for 
treating Bose-Einstein condensation in a mean-field picture, gauge symmetry breaking is 
the necessary and sufficient condition. This also concerns any perturbation theory starting 
with a mean-field approximation. 

5 Conclusion 

A Bose system with Bose-Einstein condensate is considered. The properties of the struc- 
ture factor are investigated in the frame of the self-consistent mean-field approximation, 
which is both conserving and gapless. The completed analysis demonstrates the impor- 
tance of preserving the approximation order. Going outside of the region of applicability 
of the used approximation can lead to the appearance of unphysical divergences. In the 
system with broken gauge symmetry, it is crucially important to accurately take account 
of anomalous averages. Omitting the latter is unjustified and results in the occurrence of 
fictitious instability. It is also shown that gauge symmetry breaking is necessary in any 
theory based on a mean-field picture. 
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